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ABSTRACT 


An  all-thruster  three-axis  stabilized  attitude  control 
system  has  been  designed  for  the  Naval  Postgraduate  School 
(NPS)  satellite  bus,  ORION.  The  satellite  is  a  cylinder,  19 
inches  in  diameter,  35  inches  in  length,  250  pounds  maximum 
mass  with  32  pounds  for  payloads.  ORION  will  be  ejected  from 
an  extended  Get-Away-Special  (GAS)  canister.  Launch  from  any 
GAS  can  configured  expendable  booster  or  the  space  shuttle  is 
assumed.  The  minimization  techniques  of  Pontryagin  have  been 
used  to  derive  control  laws  that  support  fuel  efficent 
operation.  A  minimum  time  cost  function  is  applied  in  the 
acquisition  phase  to  reduce  rates  to  acceptable  levels.  A 
weighted  minimum  fuel-time  cost  function  is  used  during  the 
on-station  phase.  Bang-Off-Bang  control  with  two  switching 
curves  is  employed  outside  of  a  boundary  region.  Inside  the 
boundary  region,  four-pulse  limit  cycle  control  with  time 
constants  on  the  order  of  100  seconds  is  applied.  Lifetimes 
in  the  range  of  4  to  28  months  are  obtained  utilizing  these 
techniques. 


iv 


THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed 
in  this  research  may  not  have  been  exercised  for  all  cases  of 
interest.  While  every  effort  has  been  made,  within  the  time 
available,  to  ensure  that  the  programs  are  free  of 
computational  and  logic  errors,  they  cannot  be  considered 
validated.  Any  application  of  these  programs  without 
additional  verification  is  at  the  risk  of  the  user. 
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I.  INTRODUCTION 


The  current  ORION  design  uses  spin  stabilization,  but  a 
truly  general  purpose  satellite  should  provide  attitude 
control  design  options  for  possible  payloads.  In  surveys  of 
likely  Navy  and  Air  Force  Space  Test  Program  (STP)  payloads 
during  the  past  decade,  the  majority  of  the  requests  were  for 
three-axis  stabilization  as  shown  in  Table  1-1  [l].  This 
bears  out  for  the  potential  customers  who  have  shown  an 
interest  in  the  Orion  bus.  The  intent  of  this  thesis  is  to 
investigate  a  mass  expulsion  three-axis  stabilization  design 
as  an  ORION  bus  option. 

A  three-axis  stabilized  Orion  bus  will  undoubtedly  be 
more  expensive  than  a  spin  stabilized  version.  Limited 
changes  to  the  original  design  are  desirable.  Modified 
deployable  solar  panels  will  be  necessary  to  achieve  the 
required  end-of-life  power.  Additionally,  a  revised  sensor 
package  with  more  complex  sun,  horizon  sensors  and  rate  gyros 
will  be  reviewed.  This  thesis  investigates  minimum  fuel-time 
optimization  with  respect  to  thruster  operation  using 
Pontryagin's  Minimization  Principle.  The  effects  of  external 
disturbances  including  gravity  gradient  and  aerodynamic  drag 
cannot  be  ignored.  Models  of  a  three-axis  ORION  were 
simulated. 
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Table  1-1  SUMMARY  OF  TYPICAL  NAVY/STP  PAYLOAD 
REQUIREMENTS  [1] 


Mass 

32  lbm 

Voxume 

2.36  ft3 

Power 

34  watts 

DataRate 

5000  bits/sec 

Orbit 

200-800  nm  circular 

Inclination 

0“-30*  or  60 0 -120 ° 

Instruments 

Particle  counter  or  Lens 

Control 

3  Axis, ±.75° 
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II.  BACKGROUND 


Today  most  satellites  built  in  Western  countries  are 
large,  complex  and  multi-mission/payload  capable.  The  cost 
of  additional  payloads  has  proven  exponential  due  to  many 
factors,  of  which  the  most  debilitating  is  that  contractors 
are  not  motivated  to  organize  their  program  structures 
efficiently.  Other  factors  include:  expensive  launch  from 
only  two  available  sites,  the  loss  incurred  when  a  booster  or 
satellite  subsystem  fails  and  the  increased  testing  required 
before  launch.  This  puts  the  free  world  in  a  highly 
undesirable  position.  Reserve  assets  and  the  ability  to 
produce  large  numbers  quickly  i.e.,  surge  capability  for  these 
high  value  satellites  is  extremely  limited,  if  not  non¬ 
existent.  In  addition,  most  programs  have  been  stalled  at  the 
1970-75  technology  level  due  to  a  general  fear  of  using 
unproven  components  when  the  probability  of  success  is  already 
questionable  and  the  cost  is  unacceptable.  The  concept  of 
small,  cheap  satellites  should  not  work  against  but  with 
current  programs.  A  multi-site  launch  capability  exists  from 
ships,  submarines,  aircraft  and  mobile  platforms.  New 
technology  can  be  implemented  and  space  qualified  on  low  value 
vehicles.  The  ability  to  build  up  reserves  and  to  launch 
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several  satellites  at  a  time  becomes  viable.  The  advantages 
are  myriad. 

A.  CURRENT  SMALL  SATELLITES 

The  present  U.s.  menu  of  small  proven  space  platforms 
consists  of  SPARTAN,  NUSAT ,  GLOMR  and  LDEF  [1].  SPARTAN  is 
a  retrievable  rectangular  platform  launched  from  the  Space 
Shuttle  that  can  fly  alongside  for  only  2.5  days  with  a 
payload  limit  of  1000  pounds.  NUSAT  (Northern  Utah  Satellite) 
is  also  Shuttle  deployed  but  from  a  Get-Away-Special  (GAS) 
canister  with  a  specially  designed  launch  mechanism  displayed 
in  Figure  2-1. 
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GLOMR  (Global  Message  Relay)  was  a  DARPA  (Defense  Advanced 
Research  Project  Agency)  payload  which  used  the  NUSAT  design. 
LDEF  (Long  Duration  Exposure  Facility)  is  basically  a  large 
(up  to  20,000  pounds) ,  multi-experiment  platform  that  provides 
a  facility  for  extended  space  exposure.  Figures  2-2,  2-3  and 
2-4  illustrate  SPARTAN,  NUSAT  and  GLOMR.  Though  these 
represent  advancements  in  small  satellite  development,  they 
are  fundamentally  free-flying  structures  incapable  of  transfer 
to  higher  orbits  or  long  duration  attitude  control  [1]. 


Figure  2-2  Structure  for  Spartan  Spacecraft  [1] 
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Figure  2-4  DARPA  GLOMR  Satellite  [1] 
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The  ability  to  produce  small  multi-purpose  satellites  has 
existed  since  the  dawn  of  the  space  age.  The  need  to  produce 
a  simple  yet  reliable,  low  cost  satellite  bus  that  is 
available  to  commercial  and  military  interests  alike  has  never 
been  as  great  as  it  is  now.  The  Naval  Postgraduate  School 
(NPS)  satellite  bus  concept  not  only  encompasses  propulsion 
and  attitude  control  subsystems,  but  also  standard  electric 
power,  telemetry,  thermal  and  computer  subsystems. 

B.  ORION  HISTORY 

The  history  of  NPS  satellite  development  began  in  spring 
1985  when  the  then  Chairman  of  the  NPS  Space  Systems  Academic 
Group,  Dr.  Allen  Fuhs,  taught  the  annual  spacecraft  design 
course.  The  project  for  that  particular  class  of  space 
engineering  and  operations  students  was  the  design  of  a  small 
submarine  launched  satellite  [1].  The  course  was  taught  two 
more  years  before  Dr.  Fuhs  retired  from  NPS.  Each  of  these 
two  classes  worked  on  applications  of  the  satellite.  One 
student,  LCDR  Austin  Boyd,  was  so  intrigued  with  the  basic 
concept  that  he  pursued  a  thesis  focused  on  the  simple, 
reliable  design  of  three  satellite  subsystems:  structure, 
propulsion  and  attitude  control.  LCDR  Boyd  proposed  a  Shuttle 
GAS  canister  deployable  satellite,  like  GLOMR  or  NUSAT.  The 
GAS  program  is  a  Shuttle  option  that  allows  inexpensive, 
autonomous  experiments  to  fly  space  available  in  the  cargo 
bay  [1].  Due  to  the  unceasing  efforts  of  LCDR  Boyd  and  Dr. 
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Fuhs,  agencies  such  as  DARPA,  SDIO  (Strategic  Defense 
Initiative  Organization)  and  Navy/AiT  Force  STP  (Space  Test 
Program)  Offices  began  showing  a  marked  interest  in  the  GAS 
satellite  named  ORION  and  in  small  spacecraft  in  general. 
This  interest  prompted  the  USAF  STP  office  to  pursue  the 
development  of  an  extended  and  enlarged  GAS  canister  with  an 
improved  launch  mechanism.  The  improved  GAS  can  illustrated 
in  Figures  2-5  and  2-6  can  hold  a  satellite  that  is  35  inches 
long,  19  inches  in  diameter  and  can  weigh  up  to  250  pounds 
[1].  In  addition  to  the  NPS  satellite,  three  other  efforts 
were  launched  in  1986  toward  small,  general  purpose  design. 
These  included:  "T-SAT"  by  INTRASPACE  Company,  an  advanced 
NUSAT  by  GLOBESAT  and  American  Rocket  Company  and  an  Orion 
type  design  by  TRW  [1].  This  sudden  flurry  of  activity  in  the 
small  satellite  business  was  capped  by  an  encouraging  DARPA 
program  initiated  in  early  1988  called  Lightsat.  The  Lightsat 
program  promotes  low-cost,  lightweight  satellites  to  be 
utilized  by  operational  field  commanders  [1].  It  appears  that 
the  idea  of  the  small,  general  purpose  satellite  has  finally 
come  of  age. 

C.  ORION  PRELIMINARY  DESIGN 

The  preliminary  design  of  ORION  adopted  by  LCDR  Boyd  is 
simple,  reliable  and  relatively  inexpensive  with  a  price  tag 
of  about  $1.5  million.  "A  general  purpose  architecture  was 
defined  as  the  ability  to  successfully  integrate  various 
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payloads  of  the  proper  size  while  providing  a  propulsion, 
attitude  control  and  standardized  electrical,  data  and 
attitude  control  interfaces.'*  [1] 


Figure  2-5  USAF/BALL  Aerospace  Extended  GAS  Canister  [1] 


General  purpose  naturally  implies  that  trade-offs  must  be 
made  in  the  design  of  the  payload  in  order  to  reap  the 
benefits  of  a  low  cost  bus.  Trade-offs  are  weight  and  size 
driven.  The  payload  of  any  generic  bus  is  constrained  to 
weight  and  volume  limits  that  do  not  affect  most  designs  where 
the  bus  is  built  around  the  payload.  A  brief  description  of 
ORION  subsystems  will  allow  insight  into  the  preliminary 
design. 

The  ORION  structural  subsystem  weighs  approximately  40  lbm 
and  consists  of  a  .62"  thick  baseplate,  three  .75"  honeycombed 
equipment  decks,  four  extruded  1/16"  thick  longerons  which 
house  78.5"  long  booms,  a  1"  thick  honeycombed  propellant  tank 
strongback  and  four  .05"  thick  semicircular  skin  panels.  The 
baseplate,  longerons,  booms  and  skin  panels  are  made  of  7075- 
T6  aluminum  while  the  honeycombed  items  use  stainless  steel 
facing  material.  The  booms  are  designed  to  support 
magnetometers  or  2  lbm  tip  weights  to  provide  stability  about 
the  longitudinal  spin  axis.  During  Shuttle  launch,  this 
cylindrical  structure  will  experience  less  than  .035"  of 
deflection.  In  compliance  with  NHB  1700. 7A  safety 
specifications,  the  ORION  model  resonates  above  35  Hz  at  modes 
of  160.8,  178.3  and  244.6  Hz  [1].  Figure  2-7  demonstrates  the 
basic  structural  subsystem. 

The  propulsion  subsystem  uses  hydrazine  monopropellant  and 
Shell  405  catalyst.  The  hydrazine  tank  is  a  spherical 
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Figure  2-7  ORION  Structural  Mock-up  [1] 


positive  expulsion  type.  It  is  nitrogen  pressurized  and  can 
hold  up  to  71.5  lbm  of  monopropellant.  One  40  lbf  thruster 
will  be  used  for  orbital  maneuvers  while  six  0.1  lbf  thrusters 
will  be  used  for  attitude  control  [l].  Design  requirements 
included  the  ability  to  transfer  from  a  Shuttle  orbit  of  135 
run  to  835  nm,  the  lower  limit  of  the  Van  Allen  radiation 
belts,  in  addition  to  a  lifetime  expectancy  of  90  days  to  3 
years  [1]. 

The  ORION  power  subsystem  consists  of  1765  in2  of  the 
outer  skin  covered  by  solar  cells  and  a  back-up  supply  of  24 
NiCad  cells  rated  at  1.25  volts  each  [1].  An  end-of-life  DC 
requirement  of  60  watts  is  based  on  the  following  assumptions: 
(1)  one  third  of  the  skin  surface  area  will  be  exposed  to  the 
sun,  (2)  solar  cell  efficiency  of  14%,  (3)  solar  constant  of 
125.6  watts/ ft2  and,  (4)  10%  degradation  per  year.  Even  with 
these  liberal  assumptions,  end-of-life  power  occurs  after  only 
18  months. 

The  ORION  attitude  control  subsystem  consists  of  six  0.1 
lbf  thrusters,  one  sun  sensor,  one  horizon  sensor  and  four 
magnetometers  mounted  on  78.5"  booms.  Spin  stabilization 
about  the  longitudinal  axis  with  ±1  degree  of  pointing 
accuracy  was  chosen  as  the  simplest,  most  cost  effective 
attitude  control  design  option  for  ORION  [1].  The  fuel  budget 
of  the  attitude  control  subsystem  is  limited  by  transfer 
altitude  and  effective  deployment  of  the  stabilizing  booms  in 
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this  configuration.  Figure  2-8  describes  the  attitude  control 
and  propulsion  subsystems. 
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The  remaining  ORION  subsystems  were  briefly  discussed  by 
LCDR  Boyd.  The  thermal  control  subsystem  consists  of  strip 
heaters  and  insulating  blankets.  The  health  and  welfare  of 
Orion  will  be  monitored  by  a  computer  subsystem  and  telemetry 
link  [1].  Note  that  the  preliminary  design  assumes  Space 
Shuttle  deployment  which  is  questionable  due  to  the  as  yet 
approved  use  of  hydrazine  in  a  GAS  can.  After  all  the 
subsystem  masses  are  tallied,  a  total  of  about  32  lbm  remains 
for  payload  usage.  Figure  2-9  shows  the  component  placement 
for  the  ORION  preliminary  design. 


Figure  2-9  ORION  Cross- 
section  [1] 


Table  2-1  ORION  COMPONENT  MASSES  [1] 


Component  Mass ( lbm) 


Top  plate 

1.65 

Longerons (4) 

5.80 

Payload  Midplate 

1.65 

Propellant  Tank  Strongback 

4.20 

Baseplate 

9.75 

Launch  Restraint  Pins (8) 

.80 

Pressurant  Mid-deck 

1.05 

Batteries/Boxes ( 4 ) 

28.0 

Earth  Sensor 

.  30 

Propellant  Tank 

13.0 

Hydrazine  Control  Thrusters (6) 

9.18 

Orbital  Transfer  Thruster 

1  93 

Attitude  Control  and  Payload  Computer 

16.0 

Telemetry  Transmitter/Receiver ( 2 ) 

10.0 

Pressurant  Gas  and  Bottle (2) 

5.00 

Fill  and  Drain  Valve(3) 

2.10 

Assorted  Propulsion  Tubing 

2.00 

Pyrotechnic  Valve (2) 

1.00 

Power  Conditioning  Circuits (2) 

4.00 

Payload 

32.0 

Skin  and  Fasteners 

12.0 

Solar  Cells 

.70 

Booms (4) 

4.20 

Hydrazine  Fuel 

71.5 

Magnetometers ( 4 ) 

2.00 

Total  Mass 

245.81 
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III.  SPECIAL  THREE-AXIS  CONSIDERATIONS 


A  three-axis  stabilized  ORION  will  impose  additional 
constraints  on  the  structural,  propulsion  and  power  subsystem 
design  compared  to  a  spin  stabilized  version.  A  major 
consideration  is  the  increased  fuel  usage  for  a  mass  expulsion 
attitude  control  system  (ACS) .  There  are  several  options 
available  to  extend  the  ORION  lifetime.  Although  not 
discussed  in  this  thesis,  but  promoted  as  avenues  of  further 
research,  there  exist  other  energy  efficient  three-axis 
attitude  control  schemes.  [2] 

A.  ALTERNATE  THREE-AXIS  SCHEMES 

Momentum  bias  systems  consist  of  a  momentum  wheel  and 
either  mass  expulsion  thrusters  or  magnetic  torquers  to  dump 
excess  momentum  and  control  motion  about  the  other  two  axes. 
The  momentum  wheel  supplies  gyroscopic  stiffness  which,  in 
essence,  cancels  unwanted  moments  about  one  axis.  Significant 
fuel  savings  result  in  the  thruster  moment  bias  design  in 
comparison  to  the  all-thruster  system,  but  the  added  mass  of 
the  wheel  must  be  compensated  to  maintain  the  same  payload 
mass.  Magnetic  torquers  use  the  torque  produced  when  a 
magnetic  dipole  cuts  through  the  Earth's  magnetic  field,  but 
they  are  limited  to  low  Earth  orbits.  They  provide  an  endless 
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reservoir  of  control  energy  as  long  as  current  flows  in  the 
satellite  magnets,  but  they  are  not  preferred  if  fast  reaction 
times  are  required.  The  necessity  for  propellent  is 
diminished  unless  orbital  transfer  is  desired.  [2] 

Gravity  gradient  boom  and  thrusters  or  magnetic  torquers 
afford  a  less  expensive  form  of  control  since  the  boom  is 
passive.  Accuracy  is  reduced  compared  to  the  active  momentum 
bias  system,  but  fuel  is  still  conserved. 

Reaction  wheel  systems  operate  independently  about  each 
axis  just  as  a  three-axis  mass  expulsion  system  does.  [2] 
Thrusters  are  employed  to  dump  excess  wheel  momentum  so  a  fuel 
savings  is  realized  when  equated  to  the  all-thruster  system. 
The  disadvantages  of  reaction  wheel  systems  include  the 
increased  weight,  complexity  and  cost  of  three  wheels  which 
could  disqualify  this  system  for  a  small  satellite  of  ORION 
caliber. 

B.  FUEL  SAVING  METHODS 

In  addition  to  the  minimum  fuel-time  control  laws 
discussed  and  simulated  in  Chapters  V  and  VI,  fuel 
conservation  can  be  practiced  by  restricting  the  employment 
of  the  40  lbf  orbital  insertion  thruster.  Confining  usage  of 
the  40  lbf  thruster  to  negation  of  the  effects  of  aerodynamic 
drag  could  liberate  most  of  the  hydrazine  for  on-station 
attitude  control. 
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The  ability  to  transfer  between  circular  orbits  is  a 
luxury  for  a  small  satellite,  but  AV  required  is  lessened  when 
the  booster  that  carries  ORION  can  insert  close  to  the 
requested  altitude.  The  original  ORION  designers  assumed  that 
the  Space  Shuttle  would  be  the  only  practical  launch  option. 
In  light  of  the  stringent  safety  requirements  for  a  man-rated 
payload,  the  shuttle  seems  to  be  the  least  feasible  launch 
alternative  when  dealing  with  hydrazine  monopropellant.  Table 
3-1  lists  current  launch  vehicles,  payload  weights,  altitude 

and  launch  costs  [6]. 

Table  3-1  SMALL  SATELLITE  BOOSTER  OPTIONS 


Booster  Name  Payload (kg)  Altitude (km)  Cost(106) 


Titan 

14,030 

LEO 

$90+ 

4525 

GEO 

Atlas/Centaur 

5882-7625 

LEO ( i=28 . 5  ° 

$40 

2330-3145 

GTO ( i=28 . 5 0 

Delta  II 

4460 

LEO 

$35-50 

1891 

GTO 

Scout  G-l 

219 

555 

$11-12 

*Minuteman  III 

680 

LEO (polar) 

$5 

♦Conestoga  II 

724 

LEO ( i=98  * ) 

$7.03 

1040 

LEO ( i=38  ° ) 

♦Super  Starbird 

2080 

LEO ( i=0 * ) 

$6+ 

724 

LEO (polar) 

Liberty  1A 

7260 

LEO 

$2-4 

Pegasus 

407 

LEO ( i=0  ° ) 

$6 

271 

LEO (polar) 

Foreicrn  Svstems 

Ariane  IV (France 

4200 

GTO 

$80+ 

Long  March (China) 

1267-2489 

LEO 

$35+ 

271-1358 

GTO 

♦H-2 (Japan) 

8978 

LEO 

7 

1990 

GEO 

Proton (USSR) 

20, 000 

LEO ( i=52  0  ) 

$25-40 

2000 

GEO 

*not  commercially  available  yet 
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If  ORION  must  transfer  to  a  higher  orbit,  the  fuel  budget 
will  be  driven  by  the  AV  required.  The  equation  for  the 


velocity  of  a  circular  orbit  is 


Vc  =  (^/R)%  (3.1) 

where  /i=GMe=3 . 986xl05  km3/s2,  g  is  the  Universal  gravitational 
constant  and  M,  is  Earth's  mass.  The  radius  R=R.+h,  R.=6378  km 
is  Earth's  radius  and  h  is  the  satellite  altitude.  A  Hohman 
transfer  is  the  most  efficient  means  of  transitioning  between 
coplanar  orbits  and  is  depicted  in  Figure  3-1. 


The  total  AV  required  for  a  Hohman  transfer  is 

AVi  =  Vi  -  Vcl 

av2  =  vc2  -  v2 

AVt  =  AV:  +  AV2  (3.2) 
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where  AV:  is  the  difference  between  the  elliptical  transfer 
orbit  velocity  at  perigee  and  the  inner  circular  orbit 
velocity  and  aV2  is  the  difference  between  the  transfer  orbit 
velocity  at  apogee  and  the  desired  circular  orbit  velocity. 
The  elliptical  orbit  transfer  velocity  is 

V  -  ( ( 2/i/R)  +  E)%  (3.3) 

The  specific  mechanical  energy  of  the  transfer  orbit  is 

E  =  -  p/2a  (3.4) 

where  a  is  the  length  of  the  semi-major  axis  of  the  transfer 
ellipse.  The  thrust  equation  for  a  rocket  engine  is 

F  =  meVe  +  Ae(pe  -  p0)  (3.5) 

me  =  mass  flow  rate 
Ve  =  exhaust  velocity 
Ae  =  exhaust  cross-section 
pe  =  exhaust  pressure 
p0  =  ambient  pressure 

Specific  impulse  is  a  performance  measure  of  the  thrust  to 
mass  flow  rate  ratio  for  a  rocket 

ISP  =  F/m  (3.6) 
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ISP  is  a  function  of  propellant  type.  The  specific  impulse 
of  hydrazine  is  ISP=235  sec.  The  total  impulse  is  the 
aggregate  thrust  obtained  over  time 


IT  =  J\  Fdt  =  ISP  mp  (3.7) 


where  11^  is  propellant  mass.  The  total  propellant  mass 
required  to  boost  a  satellite  to  a  higher  orbit  using  a  Hohman 
transfer  is 

mp  -  mQ  (1-  e<-iV'ISP‘>)  (3.8) 


where  m0  is  the  original  satellite  mass.  The  curves  in  Figure 
3 -2  demonstrate  the  propellant  mass  required  to  boost  ORION. 


Figure  3-2  Hydrazine  Required  by  ORION  for 
Transfer  Orbits 
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ORION  will  only  carry  71.4  lbm  of  hydrazine  so  the  fuel 
remaining  after  transfer  is  alloted  for  attitude  control.  [1] 


C.  SOLAR  PANELS 

Three-axis  stabilized  satellite  power  subsystems  generally 
include  arrays  of  solar  panels  that  extend  after  the  transfer 
orbit.  Panels  allow  more  direct  utilization  of  the  sun's 
energy  by  maximizing  the  perpendicular  rays  of  received 
sunlight  through  proper  panel  rotation.  Solar  cell  output 
voltage  is  optimized  in  this  configuration.  This  is  the  topic 
of  another  thesis  design  study  and  will  not  be  examined  in 
detail  in  this  or  later  chapters.  Solar  panels  will  add  more 
mass  to  the  ORION  structure  and  thereby  add  further  size 
constraints  to  other  subsystems.  Suggested  solar  array 
designs  for  a  three-axis  ORION  consist  of  cylindrical  body 
mounted  arrays  that  deploy  or  accordion  style  arrays  that 
unfold  and  extend  as  shown  in  Figures  3-3  and  3-4. 


Mounted  Solar  Array  [5] 


Solar  panels  alter  the  moments  of  inertia  of  a  satellite 
when  deployed.  ORION  attitude  control  design  assumes  that 
IX=IZ  which  greatly  simplifies  the  control  problem  by  causing 
one  axis  to  act  independently  of  the  other  two  in  all  regimes 
of  motion.  For  Ix*Iy*It,  coupling  occurs  between  all  axes  when 
the  angular  rate  errors  are  large.  An  axially  symmetric  array 
would  alleviate  the  coupling  in  one  axis. 

The  additional  constraints  on  fuel  and  power  budgets 
levied  by  a  three-axis  satellite  are  not  insurmountable,  but 
they  must  be  investigated  in  more  detail  in  future  ORION 
design  analyses. 
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IV.  THREE-AXIS  ATTITUDE  DETERMINATION 


Satellite  attitude  control  systems  (ACS)  consist  of 
angular  position  and  rate  sensors  or  estimators  and  control 
actuators  as  shown  in  Figure  4-1.  The  sensors  and  estimators 
feed  back  angle  and  angle  rate  information  in  the  form  of 
error  signals  to  be  processed  as  required  through  a  programmed 
control  law.  The  actuators  then  react  to  commands  from  the 
ground  or  on-board  microprocessors  to  implement  corrections 
in  attitude.  The  satellite  is  affected  by  unwanted  internal 
and  external  torques  such  as  gravity  gradient,  aerodynamic 
drag  and  thruster  misalignment.  An  ACS  conforms  to  the 
equations  of  motion  and  reacts  to  known  and  unknown 
disturbances.  (2] 


Figure  4-1  Attitude  Control  Block  Diagram  [2] 
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A.  SATELLITE  DYNAMICS 


In  modeling  the  ORION  ACS ,  rigid-body  dynamics  will  be 
assumed.  Rigid  body  dynamics  are  used  to  approximate  systems 
that  have  limited  flexibility.  The  equations  of  motion  for 
a  satellite  can  be  described  by  the  movements  of  a  point  mass 
in  space  as  shown  in  Figure  4-2. 


The  mass  P  can  experience  rotation  or  translation  with  respect 
to  the  inertial  reference  frame  XYZ.  The  position  vector  of 
P  is 

R  =  Rc  +  r  (4.1) 

where  R,  is  a  vector  from  the  XYZ  inertial  frame  to  the  xyz 
body  frame  and  r  is  the  position  vector  of  P  in  the  xyz  frame. 

r=  xi  +  yj  +  zk  (4.2) 
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The  absolute  velocity  of  P  is 


V  = 


R 


•  • 
Re  +  r 


(4.3) 


where  r  includes  linear  and  angular  velocity  in  the  body 
frame. 

•  • 

r  =  rrel  +  u>  x  r  (4.4) 

Substituting  (4.4)  into  (4.3),  the  absolute  velocity  becomes 

V  =  Vc  +  Vrel  +  w  x  r  (4.5) 

The  absolute  acceleration  of  P  is  the  derivative  of  absolute 
velocity. 


a=  ac  +  ar8l  +  2w  x  Vrel  +  w  x  (w  x  r)  (4.6) 

In  order  to  simplify  future  equations,  two  assumptions 
will  be  made.  The  xyz  frame  is  internal  to  the  spacecraft 
body  so  that  vt.i=o  and  the  spacecraft  center  of  mass  (CM) 
coincides  with  the  center  of  the  body  frame  so  that  /mrdm=0. 
The  linear  momentum  of  a  rigid-body  is  defined  as  the  integral 
of  the  absolute  velocity  times  the  differential  mass. 
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p  =  Xni(vc  +  u  x  r)dm  =  mV, 


(4.7) 


Angular  momentum  of  a  rigid-body  is  defined  as  the  moment  of 
the  linear  momentum  or  the  cross  product  of  r  and  P. 

Hc  =  /^(r  x  V) dm  =  /„  r  x  («  x  r)dm  (4.8) 

The  angular  velocity  vector  is  equal  to 

o>  —  wxi  +  wyj  +  wzk  (4.9) 

Moments  and  products  of  inertia  are  analogous  to  the  mass 
of  a  linear  system  but  relate  to  mass  distribution  for  a 
rotating  rigid  body.  Inertia  is  the  sum  of  the  products  of 
the  component  masses  and  their  distance  form  the  chosen  axis 
of  rotation. 

I  =  Z  m;  r i2  (4.10) 

Moments  of  inertia  become 

Ix  =  S  mjy2  +  z“)  (4.11) 

Iy  =  S  iMxz  +  z2) 

Iz  =  S  iMx2  +  y2) 
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Table  4-1  lists  preliminary  subsystem  masses  and  principal 
moments  of  inertia  for  ORION. 

Table  4-1  SUBSYSTEM  MOMENTS  OF  INERTIA  [12] 


Subsystem  Mass (kg)  lx  Iy  Iz 


Structure 

28.99 

2.4652 

2.4652 

1.4148 

Propulsion 

44.57 

1.6218 

1.6647 

0.8844 

Telemetry 

4.61 

0.6089 

0.5500 

0.0652 

Attitude  Control 

4.41 

9.3621 

9.3827 

18.6291 

Computer 

3.45 

0.3482 

0.3482 

0.0214 

Electrical  Power 

13.84 

0.7379 

0.5016 

0.3011 

Thermal 

0.0 

0.0000 

0.0000 

0.0000 

Payload 

11.22 

1.0812 

1.1117 

0.0916 

Total  (kg-m2) 

111.09 

16.23 

16.02 

21.04 

Total  (Slug-ft2) 

7 . 61 

11.9 

11.82 

15.97 

Without  booms 
Total  (kg-m2) 

106.68 

6.8633 

6.6415 

2.7758 

Total  (Slug-ft2) 

7.31 

5.06 

4.90 

2.05 

The  products  of  inertia  below  are  zero  when  the  principal 
moments  of  inertia  are  chosen  as  rotation  axes. 


I*y  =  s  m,  xy 
Ixz  =  Z  mA  xz 
Iyz  =  S  mj  yz 

After  solving  (4.8)  and  recognizing  (4.9)  and  (4.11),  the 
angular  momentum  becomes 


ix  ~ ixy  "”ixz 

wx 

Hc  = 

~ixy  iy  “  iyz 

_  —  ixz  ~ iyz  iz  . 

.wz. 
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The  inertia  matrix  is  called  the  inertia  tensor.  To  simplify 
the  design  of  ORION,  the  inertia  tensor  will  be  assumed 
diagonal  with  zero  cross  terms. 

The  kinetic  energy  of  a  rigid  body  is  defined  as  the 
integral  of  the  dot  product  of  the  absolute  velocity. 

T  =  i/B  V-V  dm  (4.13) 

This  thesis  will  focus  on  rotation  only.  Substituting  (4.5) 
into  (4.13)  the  kinetic  energy  reduces  to 


“i  wy  W2  ] 

H 

X 

O 

O 

_ i 

o 

>> 

H 

O 

Oy 

O 

O 

H 

N 

.wz. 

(4.14) 


Newton's  Second  Law  can  be  applied  to  a  rigid  body  as 

F  =  /„  adm  =  mac  (4.15) 

A  moment  about  the  satellite  CM  is  the  cross  product  of 
distance  and  force. 


Mc  =  Hc  +  /mRc  x  rdm  =  Hc 

=  Hcrel  +  «  x  Hc 


(4.16) 


Angular  momentum  in  the  body  frame  can  be  expressed  as 
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Hc=  Hxi  +  Hyj  +  Hzk 


(4.17) 


where  xyz  coincides  with  the  principal  moments  of  inertia. 
From  (4.12),  the  elements  of  the  angular  momentum  vector  are 

H*  =  (4.18) 

Hy  =  IyWy 
Hjj  ^2^2 


After  some  manipulation,  (4.16)  is  transformed  into  the  Euler 
Moment  equations. 


Mx  = 

X 

3 

K 

H 

+ 

wywz  (  ^*z 

-  Iy) 

My  = 

• 

Iywy 

+ 

WxWzUx 

-  I.) 

Mz  = 

• 

^2^2 

+ 

WxWy(Iy 

-  I*) 

(4.19) 


Equation  (4.19)  will  be  used  extensively  through  out  this 
thesis  to  simulate  ORION  attitude  control. 

An  inertial  coordinate  system  with  X0Y0Z0  fixed  at  the 
center  of  the  Earth  as  shown  in  Figure  4-4  is  used  as  a 
reference  frame  when  describing  satellite  motion.  The  roll 
axis,  X,  describes  the  orbital  velocity  vector  with  magnitude, 
wQ,  which  is  assumed  constant  for  a  nearly  circular  orbit. 
The  yaw  axis,  Z,  extends  from  the  satellite  CM  toward  the 
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Figure  4-4  Spacecraft  Coordinate  System  [2] 

Earth  CM.  The  pitch  axis,  Y,  is  perpendicular  to  the  orbital 
or  XY  plane.  Errors  caused  by  disturbance  torques  are 
described  by  the  angles  <f> ,  j/>  and  6  about  the  roll,  yaw  and 
pitch  axes,  respectively.  In  unperturbed  form,  the  Euler 
moment  equations  become 


The  right  hand  side  of  (4.20)  exhibits  the  coupling  that 
occurs  when  angular  error  rates  are  large.  For  on  orbit 
control  of  a  three-axis  stabilized  system  the  angular 
velocities  are  assumed  small  which  reduces  (4.20)  to 

(4.21) 


Satellite  acquisition  is  a  stage  in  which  initial  large 
attitude  adjustments  are  made  to  correct  for  booster 
separation.  In  this  phase  the  satellite  frequently 

experiences  a  tumbling  motion  which  translates  to  large  errors 
in  angular  position  and  rate.  Excess  rates  must  be  negated 
before  the  on  orbit  control  law  can  be  implemented.  Angular 
velocity  errors  as  viewed  from  the  inertial  frame  can  be 
expressed  as 


Ixwx  =  0 


:LyWy  =  o 


Izu>z  =  0 


wx  =  <t>  ~  w0^  (4.22) 

wy  =  6  -  w0 

-  ij)  +  u0<}> 

Replacing  rate  errors  in  the  Euler  moment  equations  yields 
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•  •  •  • 

M*  =  I M  -  +  (I,  -  +  <*o4>)  (4.23) 

My  =  Iy  9 

Mz  =  Iz(tf>  +  W0^)  +  (Iy  -  Ix)w0(^  -  w0\i>) 

B.  EXTERNAL  DISTURBANCES 

The  major  external  disturbances  for  low  altitude  orbits 
are  caused  by  the  Earth's  gravitational  field  and  aerodynamic 
drag.  Solar  pressure  contributions  are  considerd  negligible 
at  low  altitude  but  a  dominant  factor  beyond  geosynchronous 
altitudes. 

1.  Gravity  Gradient  Torque 

Gravity  gradient  torques  can  be  used  as  a  passive  form 
of  three-axis  stabilization  for  Earth  pointing  payloads  if  a 
boom  is  attached  to  the  satellite.  This  torque  results  from 
the  interaction  between  the  earth-satellite  masses.  [2]  The 
gravitational  moment  can  be  expressed  as 

He  =  J,  r  x  F  dm  (4.24) 

where  the  gravitational  force  is 

F  =  ftllk  ~  r>dBlf1  +  2E3  (4.25) 

Re3  Re"  . 

Substituting  (4.24)  into  (4.23)  results  in 
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(4.26) 


Mg  =  ~  J(r  x  Re)  (r  •  R.) 
Re' 


where  R.,  the  vector  between  Earth-satellite  centers  of  mass 
as  shown  in  Figure  4-5  is  expressed  in  the  satellite  body 
frame  as 


R„  =  R„(-sinfli  +  sin^cosSj  +  cos^cosflk)  (4.27) 
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For  small  angle  corrections  equation  (4.26)  reduces  to 


Mg  « 


*(I,  - 

•  6  (Iz  - 

0 


(4.28) 


2.  Aerodynamic  Drag  Torque 

Aerodynamic  drag  is  caused  by  the  friction  between  the 
satellite  as  it  orbits  at  an  angular  rate  of  w0  and  the 
Earth's  atmosphere.  It  is  felt  most  below  500  km.  The 
effects  of  drag  cause  a  general  decrease  in  angular  velocity 
thereby  resulting  in  decreasing  altitude  or  orbital  decay. 
The  aerodynamic  drag  force  is  a  function  of  air  density,  p; 
velocity,  V;  surface  area  perpendicular  to  the  flight  path, 
A;  and  a  dimensionless  number  called  the  coefficient  of  drag, 
CD.  [8] 


Fd  =  i  pV2CDA  (4.29) 

The  orbital  lifetime  of  a  satellite  is  a  function  of  altitude, 
h,  and  eccentricity,  e.  [7] 


L  =  NP 


1.15185X10"7  N 


6378  +  hi 


1-e 


3/2 


(4.30) 


where  N  is  the  estimated  number  of  revolutions  and  P  is  the 
period  of  the  orbit.  The  values  of  N  can  be  obtained  directly 
from  Figure  4-6  where  the  ballistic  coefficient  for  ORION  is 
m/ (CpA)  =132 . 87  kg/m2.  Note  that  the  eccentricity  for  a 
circular  orbit  is  zero,  e=0.  For  ORION,  orbital  lifetime  is 
shown  in  Figure  4-7  as  a  function  of  altitude.  The  linear 
velocity  of  a  satellite  equals 

v  =  0>0R  =  (m/R)%  (4.31) 


Figure  4-6  Orbital  Decay  [7] 


Aerodynamic  Drag 

where  R  =  R.  +  h,  the  altitude  above  Earth  CM.  The  drag 
coefficient  is  between  1  and  2  in  SI  units.  The  value  used 
for  ORION  will  be  CD=2 .  The  surface  area  of  ORION 
perpendicular  to  the  flight  path  is  A=.429m2.  This  reduces 
(4.29)  to 


Fd  =  1.71  x  lO^p/R  (4.32) 

ORION  drag  torque  will  can  be  summed  up  as 

Tdrag  =  FD  L  costf  (4.33) 
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where  L=.076m  is  the  ORION  center  of  pressure  location  along 
the  y  axis. [8]  The  effects  of  aerodynamic  drag  and 
gravitational  gradient  are  simulated  for  a  three-axis  ORION 
in  Chapter  V. 

C.  SENSORS  AND  ACTUATORS 

The  ACS  sensor  suite  is  dictated  by  the  type  of 
stabilization,  pointing  accuracy,  altitude,  orbital  parameters 
and  cost  constraints  of  the  satellite  budget.  There  are 
potentially  three  distinct  stages  of  attitude  control  that 
ORION  must  handle.  These  include  transfer  orbit,  acquisition 
and  small  error  correction  which  includes  altitude 
maintenance.  If  a  transfer  orbit  is  required,  the  satellite 
will  most  likely  be  spinning  during  this  phase  so  different 
types  of  attitude  determination  sensors  will  be  needed  to 
accommodate  spinning  and  non-spinning  modes.  In  the 
acquisition  stage,  the  satellite  experiences  large  error  rates 
which  must  be  negated  before  accurate  pointing  can  be 
sustained.  The  small  error  stage  is  the  steady  state 
operation  that  the  satellite  will  undergo  for  the  majority  of 
its  active  life. 

The  attitude  of  a  spacecraft  in  inertial  space  can  be 
completely  realized  by  two  vectors  as  long  as  they  are  not 
col inear.  Point  sources  such  as  the  sun  and  stars  afford 
excellent  reference  vectors  while  the  Earth  and  planets 
present  extended  sources.  Pitch,  roll  and  yaw  attitude  errors 


relative  to  the  desired  satellite  orientation  must  be 
determined  in  addition  to  error  rates  about  these  three  axes. 
Combinations  of  horizon  and  sun  sensors  can  furnish  position 
error  while  rate  sensing  gyros  or  error  estimation  routines 
can  be  used  to  resolve  rate  errors. 

The  current  selection  of  attitude  sensors  for  a  spinning 
ORION  consists  of  sun  sensors,  horizon  sensors  and 
magnetometers.  The  sun  and  horizon  sensors  provide  the 
primary  angular  error  output  while  the  magnetometers  act  as 
back-up  sensors.  The  spinning  motion  of  the  original 
satellite  simplifies  the  number  and  cost  of  primary  sensors 
by  producing  a  periodic  search  motion.  In  a  three-axis 
stabilized  system  this  scanning  motion  must  be  integrated  into 
the  instrument  in  the  case  of  the  horizon  sensor  while 
multiple  sun  sensors  with  wider  fields  of  view  must  be 
employed. 

1.  Horizon  Sensors 

Earth  or  horizon  sensors  are  infrared  devices  that 
detect  the  14-16  /im  radiance  profile  of  the  earth  as  the  space 
to  earth  and  earth  to  space  transitions  occur.  Horizon 
sensors  can  be  used  to  ascertain  pitch  and  roll  errors.  [2] 
Scanning  or  radiation  balance  detectors  are  employed  in  three- 
axis  designs  to  provide  a  sweeping  action  across  the  disk  of 
the  Earth.  Scanning  sensors  use  two  beams  centered  on  a 


42 


reference  line  to  scan  the  Northern  and  Southern  hemispheres 
in  an  east-west  direction.  The  scan  pulse  width  is 


pw=  d/(2hwh)  (4.34) 

where  d  is  the  distance  between  horizon  crossings,  h  is  the 
satellite  altitude  and  wh  is  the  scanning  rate  of  the  sensor. 
The  difference  between  in-crossing  to  reference  and  reference 
to  out-crossing  pulse  widths  for  each  beam  can  be  used  to 
determine  pitch  error  while  the  difference  between  Northern 
and  Southern  scans  is  proportional  to  roll  error.  An 
incremental  pulse  generator  outputs  a  pulse  at  fractions  of 
a  degree  for  the  known  scan  rate.  Pulse  width  is  computed  by 
adding  pulses  from  acquisition  of  signal  (AOS)  to  reference 
then  from  reference  to  loss  of  signal  (LOS)  .  (7]  Typical 
output  from  a  scanning  horizon  sensor  is  depicted  in  Figure 
4-8. 


Figure  4-8  Output  from  a 
Scanning  Horizon  Sensor  [7] 
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2 


Sun  Sensors 


Since  two  non-coincidental  directions  are  needed  to 
measure  alignment  uniquely,  it  is  necessary  to  find  another 
reference  vector  besides  that  from  the  Earth.  The  sun  or 
another  star  are  logical  alternatives  since  they  can  be 
approximated  as  point  sources.  Stars  are  the  best  inertial 
reference  points  and  are  not  coplanar  unless  the  inclination 
is  zero  and  the  orbital  path  is  in  the  equatorial  plane.  Star 
sensors  are  orders  of  magnitude  more  complex,  expensive  and 
bulky  than  sun  sensors  which  owe  their  simplicity  to  the  close 
proximity  of  the  most  radiant  body  in  the  solar  system. 

Sun  sensors  are  the  most  commonly  used  attitude  sensor 
due  to  their  relatively  low  cost,  lightweight  and  broad 
application.  [2]  They  can  render  accurate  position  error 
updates  in  harmony  with  horizon  sensors  when  the  satellite  to 
Earth  and  satellite  to  sun  vectors  are  non-parallel.  Vector 
alignment  can  occur  up  to  twice  an  orbit  at  local  noon  and 
midnight.  Sun  sensors  are  basically  silicon  solar  cells  with 
specialized  masks  for  varying  functions.  Analog  and  digital 
sensors  are  available  to  measure  the  angle  between  the  sun 
line  and  the  pitch  axis  or  sun  presence  in  a  particular  field 
of  view  for  instrument  protection.  The  sun  line  vector  is 
defined  as 

S  =  sinacosfil  +s.in6J  cosaCosSK  (4.35) 
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where  6  is  the  sun  declination  and  a  is  the  orbital  angle  of 
the  satellite  measured  form  local  noon.  Yaw  errors  about  the 
pitch  axis  in  the  inertial  frame  can  be  expressed  using 

j  =  -^I  +  J  +  <f>K  (4.36) 

for  small  angles.  The  sun  angle,  &a,  is  obtained  by  taking 
the  dot  product  of  the  sun  vector  and  j 

cos0s  =  S  ■  j  =  -V>sinacos5  +  sins  +  ^cosacosS  (4.37) 

As  the  yaw  error  changes  the  sun  angle  varies 


a  cos  es 

a  v> 


=  -sinacos 6 


(4.38) 


thus  yaw  error  measurement  is  obtained.  [2] 

3 .  Magnetometers 

Magnetometers  are  not  a  required  sensor  for  current 
or  future  ORION  designs.  The  original  design  utilizes  booms 
with  2  lbm  weights  at  the  end  to  enhance  the  stability  of  the 
spinning  satellite  by  altering  the  principal  moments  of 
inertia.  Three  of  the  tip  weights  are  magnetometers  which 
will  provide  a  back-up  attitude  determination  capability  that 
dead  weights  cannot.  Magnetometers  are  solenoid  coils  that 
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obey  Faraday's  Laws  of  Magnetic  Induction.  When  placed  in  a 
time-varying  magnetic  field,  solenoid  coils  will  produce  a 
voltage 


V  =  /  E  •  dl  =  -ANM(dB,/dt)  (4.39) 

where  E  is  the  electromotive  force,  A  is  the  cross-sectional 
area  of  the  solenoid  core,  N  is  the  magnetic  permeability  of 
the  core  and  B±  is  the  magnetic  field  component  along  the 
solenoid  axis.  The  voltage  is  proportional  to  the  magnetic 
field  magnitude  and  direction.  On-board  control  using 
magnetometers  requires  knowledge  of  the  Earth's  magnetic  field 
at  the  orbit  altitude.  Knowledge  of  the  magnetic  field  is  not 
complete  at  low  altitudes  and  too  weak  above  1000  km  since  the 
field  strength  decreases  as  the  cube  of  the  distance  from  the 
Earth's  CM.  [7] 

4 .  Rate  Gyroscopes 

In  the  acquisition  phase,  some  means  of  determining 
error  rates  is  required  before  they  can  be  cancelled.  Rate 
sensing  gyroscopes,  Luenberger  observers  or  Kalman  filters  for 
the  sun  and  horizon  sensor  output  are  viable  solutions.  Rate 
integrating  gyros  (RIG)  in  triplets  are  the  most  accurate 
instruments  used  to  detect  error  rates  and  positions  about 
three  axes.  The  output  of  a  RIG  is  the  current  required  to 
maintain  the  gyro  gimbals  at  a  null  or  zero  position.  The 
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current  is  proportional  to  the  torque  needed  to  zero  the 
rates.  The  total  angular  momentum  of  a  gyro  is 

H  =  LS  +  Io0O  (4.40) 

where  L  is  the  magnitude  of  the  angular  momentum  of  the  rotor 
about  the  spin  axis  and  Io0  is  the  magnitude  of  the  angular 
momentum  of  the  gimbal  system  about  the  output  axis.  The 
nulling  torque  is 

Tg  =  (dH/dt)  gyro  +  a)  x  H  (4.41) 

Angular  velocity  in  the  gyro  inertial  frame  is 

w  =  Wjl  +  w00  +  wsS  (4.42) 

where  I  is  the  gyro  input  axis,  S  is  the  rotor  spin  axis  and 
O  is  the  output  axis.  Figure  4-9  depicts  gyroscope  operation 
and  the  physical  location  of  the  axes  described  by  equations 
(4.40)  and  (4.42).  The  differential  equation  of  motion  for 
a  gyroscope  is 

I J  +  D6  +  K6  -  wjL  =  0  (4.43) 

where  D  is  a  viscous  damping  term  and  K  is  a  restoring  term. 
[7] 
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5 .  Hydrazine  Thrusters 

The  functions  of  station  keeping  and  attitude  error 
correction  are  performed  solely  by  control  actuators  called 
thrusters  or  reaction  jets  in  this  design.  The  thrusters  on 
ORION  use  the  monopropellant  hydrazine  with  the  Shell  405 
catalyst  for  these  functions.  The  hydrazine  storage  tank  is 
pressurized  with  nitrogen.  A  basic  hydrazine  reaction  jet  is 
pictured  in  Figure  4-10. 

The  lifetime  of  a  three-axis  stabilized  satellite  is 
more  acutely  limited  by  the  operational  usefulness  of  the 
Shell  405  catalyst  than  in  a  spin  stabilized  vehicle  where 
fuel  usage  can  be  more  accurately  predicted.  The  reaction 
jets  on  a  three-axis  stabilized  satellite  will  be  subjected 
to  intermittant  firing  or  low  duty  cycle  firing.  Thruster 


48 


Figure  4-10  Basic  Thruster  Design  [2] 


performance  is  largely  determined  by  the  temperature  of  the 
catalyst  bed.  During  steady  state  operation,  intermittent 
firing  of  the  jets  will  not  allow  the  catalyst  bed  to  maintain 
sufficient  warmth  to  facilitate  the  most  efficient  thrust 
profile  between  firings.  Typical  values  of  specific  impulse 
for  hydrazine  during  transfer  and  steady  state  operation  are 
235  and  135  seconds,  respectively.  The  specific  impulse  is 
seriously  degraded  in  the  steady  state  regime.  Ultimately, 
this  leads  to  the  conclusion  that  the  catalyst  bed  in  a  low 
duty  cycle  environment  must  be  heated  which  requires 
considerably  more  power  than  just  that  used  for  opening  and 
closing  valves  and  heating  the  hydrazine  to  prevent  freezing. 
[2] 

In  addition  to  the  catalyst  bed  temperature,  thruster 
performance  is  affected  by  impurities  in  the  hydrazine  that 
accumulate  in  the  catayst  bed  after  each  firing  and  variation 
in  the  nitrogen  pressure  as  the  propellant  is  depleted.  [ 2 ] 
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Figure  4-11  shows  the  linear  relationship  between  thrust 
obtainable  and  pressure. 


100  150 

Tank  pressure  (mN/m1  abs) 


Figure  4-11  Thrust  versus  Tank 
Pressure  [2] 


A  typical  thruster  profile  is  illustrated  in  figure 
4-12.  This  profile  assumes  that  physical  conditions  are 
acceptable,  but  that  there  is  a  finite  delay  between  the  time 
that  the  thruster  is  commanded  on  and  the  time  that  it  reaches 
maximum  thrust  and  likewise  for  the  command  off  sequence.  The 
time  from  tj  to  t2  and  t4  to  t5  is  a  function  of  the  valve 
opening  and  closing  action.  In  Chapter  V,  the  thruster 
profile  is  modeled  as  a  pulse  function  for  simplicity. 

The  current  thruster  selection  for  ORION  includes  a 
40  lbf  jet  for  orbital  insertion  and  six  .1  lbf  attitude 
control  jets.  Figure  4-13  lists  design  specifications  for 


these  thrusters. 
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K 


2  0* 
Dia 


87 

t. 


o 

Propellant 

Hydrazine 

□ 

Catalyst 

Shell  405/LCH  202 

□ 

Thrust,  Steady  Stale  (Ibf)  . 

40—  11 

□ 

Feed  Pressure  (psia) 

450—  100 

□ 

Chamber  Pressure  (psia) 

223  —  61 

□ 

Expansion  Ratio  . 

20:1 

□ 

Flow  Rate  (Ibm/sec) 

0  17  —  0  054 

□ 

Valve 

Wright  Components  Single  Seat 

□ 

Valve  Power 

50  Walts  @  28  vdc  and  75°F 

□ 

Weight  (Ibm) 

1.95 

Engine  . 

1,35 

Valve 

0  60 

Demonstrated  Performance 

□ 

Specific  Impulse  (Ibf-sec/lbm) 

235  —  203 

□ 

Total  Impulse  (Ibl-sec) 

4.350 

□ 

Total  pulses  . 

2.258 

□ 

Minimum  Impulse  Bit  (Ibl-sec) . 

0.7 1  @  450  psia  &  20  ms  ON 

□ 

Steady-Stale  Firing  (sec) 

97 

Figure  4-13a  Rocket  Research 
Model  107B  [1] 

40  lbf  Thruster 
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Design  Characteristics 

□  Propellant . 

□  Catalyst 

□  Thrust,  Steady  Sta  te  (Ibf) 

□  Feed  Pressure  (psia)  . 

□  Chamber  Pressure  (psia) 

□  Expansion  Ratio 

□  Flow  Rale  (Ibm/sec) 

□  Valve  . 

□  Valve  Power 

□  Weight  (Ibm) 

Engine . 

Valve  . 


.  .  Hydrazine 

. Shell  405 

. 0.252—0.042 

420—70 

. 370—60 

100:1 

0.001—0.0002 
Wright  Components  Dual  Seat 
.  9  Watts  Max.  at  28  vdc  and  45  'F 
073 
0.28 
0.45 


Demonstrated  Performance 

□  Specific  Impulse  (Ibt-sec/lbm) . 

□  Tolaf  Impulse  (ibf-sec) . 

□  Total  Pulses  . 

□  Minimum  Impulse  Bit  (Ibl-sec) 

□  Steady-State  Firing  (sec) 


. 227—206 

. 35.625 

410.000 

0  .001  @  150  psia  &  8  ms  ON 
. 64.800 


Figure  4-13b  Rocket  Research  .1  lbf  Thruster  - 
Model  103C  [1] 
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V.  MINIMUM  FUEL-TIME  PROBLEM 


Three-axis  stabilization  is  inherently  more  expensive  than 
spin  stabilization  due  to  the  increase  in  complexity  of 
attitude  sensors  and  the  energy  necessary  to  maintain  the 
desired  pointing  accuracy.  For  Orion  this  translates  to  more 
fuel  needed  once  on-station  than  the  original  spinning 
version.  This  limits  the  total  AV  available  for  orbit 
transfer  and  the  active  on-station  control  life  of  the 
satellite.  Therefore,  the  effective  operational  lifetime  of 
the  payload  could  be  dictated  by  fuel  constraints  before  other 
considerations.  In  order  to  conserve  fuel  and  provide  maximum 
operational  lifespan,  a  minimum  control  effort  solution  using 
Pontryagin's  minimum  principle  is  explored  in  this  section. 
[9] 


A.  ON-STATION  CONTROL  LAW 

For  small  angular  rate  and  pointing  errors,  examination 
of  one  axis  applies  to  all  three  axes.  Thrusters  are  the  sole 
control  device  for  this  design  study  and  act  independently  in 
this  case.  The  simplified  differential  equation  for  one  axis 
is 

I«  =  Tc  (5.1) 
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where  Tc  is  the  thruster  control  torque  and  I  is  the  inertia 
tensor  with  the  products  of  inertia  equal  to  zero.  The 
normalized  state  equations  for  one  axis  are 


Xj(t) 

0  1 

Xi(t)' 

O' 

* 

= 

+ 

x2(t). 

.0  0. 

.x2(t). 

1 

u(t) 


(5.2) 


where  xx  is  angular  position,  x2  is  angular  velocity  and  u  is 
the  control  torque  provided  by  the  thrusters.  The  minimum  fuel 
performance  index  or  cost  function  is 

J(u)  =  ft  (1  +A|u|)dt  (5.3) 

where  A  is  a  weighting  factor  applied  to  the  control.  As  A-*0, 
the  performance  index  becomes  the  minimum  time  solution  and 
as  A-*®,  the  performance  index  approaches  the  pure  minimum  fuel 
solution.  Since  any  minimum  effort  problem  is  a  trade  off 
between  elapsed  time  and  consumed  fuel  both  should  be 
considered  and  are  represented  by  (5.3).  The  performance 
measure  must  be  minimized  to  obtain  the  optimal  solution  for 
given  initial  conditions  and  final  states,  x(tf)=0.  Every 
value  of  the  weighting  factor  provides  a  minimal  trajectory. 
The  selection  of  the  weighting  factor  is  determined  by 
reasonable  transient  time  and  fuel  expenditure.  The 
normalized  admissible  controls  for  thruster  operation  are 
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I U  (t)  1  <  Tc/I  =  N 


(5.4) 


where  Tc  is  the  maximum  torque  available  from  the  reaction 
jets.  The  cost  function  for  thruster  operation  is  called 
•Bang-Of f-Bang'  since  the  satellite  drifts  before  a  final 
control  torque  is  applied  to  bring  it  to  the  desired  final 
attitude.  The  Hamiltonian  for  (5.3)  is  defined  as 

H(x,u,p,t)  =  1+X  | u |  +PiX2+p2u  (5.5) 

where  p(t)  are  Lagrange  multipliers.  The  costate  equations  are 

p^t)  =  -  au/ax1  =  o  (5.6) 

p2(t)  =  -  3H/3XZ  =  -px(t) 

which  have  solutions 

Pi(t)  =  c:  (5.7) 

p2(t)  =  -Cit  +  c2 

The  Hamiltonian  must  be  minimized  on  an  optimal  trajectory, 
so  the  form  of  optimal  control  is 
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U(t) 


P2  <  _ ^ 

-A  <  p,  <  A 
p2  >  A 


(5.8) 


N 

0 

-N 


where  switching  occurs  at  |p2|=A.  The  effect  that  any  known 
or  unknown  torque  has  on  the  satellite  can  be  best  illustrated 
in  the  phase  plane.  This  is  a  plot  of  angular  position  versus 
velocity  as  a  function  of  applied  torque.  There  are 
essentially  two  switching  curves  for  the  given  minimum 
trajectory.  These  correspond  to  a  curve  that  is  dependent  on 
the  weighting  factor,  A,  and  the  minimum  time  curve  as  shown 
in  the  phase  plane  representation  in  Figure  5-1. 


Figure  5-1  Double  Switching  Curve  Control 
Law 

If  the  initial  conditions  place  the  satellite  in  region  1,  a 
negative  torque  is  applied  until  tx  when  region  3  is  entered 
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and  control  is  turned  off.  Positive  control  is  applied  at  t2, 
when  the  satelite  reaches  region  2,  which  is  the  switching 
curve  for  the  minimum  time  solution.  The  reverse  is  true  if 
the  initial  conditions  place  the  satellite  in  region  2.  At 
tx  the  control  goes  off,  u(tx)=0,  and  (5.6)  becomes 

P2(tx)  —  x  —  -cxtx  +  c2  (5.9) 

The  position  continues  to  change  as  a  function  of  velocity  at 
turn-off.  After  t2,  control  is  turned  on,  u(t2)=0,  and  the 
satellite  ideally  moves  to  the  desired  position.  At  t2  (5.9) 
becomes 

p2(t2)  =  -X  —  ~cxt2  +  c2  (5.10) 

Subtracting  (5.10)  from  (5.9) 


(t2  -  tx)  =  2A/cx  (5.11) 

Assuming  that  velocity  is  constant  at  turn-off,  x2=0 

x2(tx)  =  x2(t2)  =  *2(t) 

Solving  the  Hamiltonian  at  the  switching  times  and 
substituting  (5.6)  for  px  yields 
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H(t1)  =  H  (t2)  =  0 

H(tx)  =  1  -  C1X2(t1)=0 
1 

Ci  =  - 

X2(tx) 

Substituting  (5.11)  into  (5.12) 

(t2  -  tx)  =  2XX2(tj) 


(5.12) 


(5.13) 


Solving  the  state  equation,  x1(t)=x2(t) 


Xi(t2)  -  Xi(tx)  =  x2(tx)  [tj-tJ 


(5.14) 


From  t2  to  tf  the  satellite  follows  a  minimum  time  trajectory 
where 


xx(t) 


i  xzZ(t> 

I N  | 


(5.15) 


Substituting 


(5.13) 

xx(t) 


and  (5.15) 

,  (4A  -f  1) 
2  I  N  j 


into  (5.14) 

x22(t)  =  0 


(5.16) 


which  is  the  control  law  for  this  minimum  fuel  problem. 
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Equation  (5.16)  becomes  usable  in  all  quadrants  when 
x22=x2|x2|  which  takes  the  sign  of  x2  into  account.  Figure  5- 
2  demonstrates  minimum  trajectories  for  various  values  of  A. 
The  initial  conditions  x1(0)=l,x2(0)=0  were  used  to  test  for 
minimization  of  the  selected  cost  function.  As  shown  in 
figure  5-3,  the  minimum  value  of  the  chosen  cost  function 
occurs  at  A=2.25  when  this  value  is  substituted  into  equation 


Figure  5-2  Parabolic  Switching  Curves  for 
A=0, 1, 10, 100 


Figure  5-4a  is  a  time  response  of  angular  position  and  rate. 
Figure  5-4b  is  the  time  response  of  control  torque.  Notice 
that  thruster  chatter  occurs  at  the  end  of  the  cycle,  since 
fine  control  near  the  desired  attitude  is  not  implemented. 
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Figure  5-4b  shows  the  "bang-deadzone-bang"  behavior  of  the 
control  law.  Figure  5-5  is  a  phase  plane  representation  of 
the  same  cycle. 


Figure  5-5  Phase  Plane  for  x-2.25 


B.  LIMIT  CYCLE  CONTROL 

Ideally,  the  last  thrust  period  will  take  the  satellite 
to  the  origin,  but  this  is  a  function  of  the  integration  step 
size.  A  similiar  analogy  exists  for  the  thruster  in  the  form 
of  the  minimum  impulse  bit  available.  Physically  angular 
rates  are  never  zeroed,  but  they  can  be  maintained  in  a 
region.  Satellite  motion  is  then  dictated  by  a  finite  limit 
cycle  even  in  the  absence  of  external  disturbances. 

While  in  the  limit  cycle  attitude  errors  are  reduced  below 
specified  values.  The  fine  tuning  region  is  determined  by 
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pointing  accuracy  requirements,  but  it  is  also  a  function  of 
minimum  jet  impulse  for  angular  velocity  and  the  precision  or 
noise  contribution  of  the  attitude  sensor  instrumentation. 
Minimum  jet  impulse  is  a  measure  of  the  smallest  thrust  unit 
available  when  the  control  valves  are  immediately  opened  and 
closed,  i.e.,  the  shortest  single  jet  pulse.  In  order  to 
increase  the  coasting  period  in  the  limit  deadzone,  the 
satellite  rate  must  be  reduced  to  the  lowest  value  possible 
without  causing  thruster  chatter.  Chatter  is  a  constant 
cycling  of  opposing  thrusters  when  limits  are  exceeded  after 
each  firing.  The  limit  cycle  for  the  cost  function  simulated 
in  this  example  is  depicted  in  Figure  5-6. 


Figure  5-6  Limit  Cycle 
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A  third  switching  curve  was  implemented  in  the  boundary 
region  only  to  minimize  fuel  expended  while  cycling.  The 
third  switching  curve  has  a  flatter  slope,  A=250,  and  thus  a 
wider  deadzone.  The  best  deadzone  obtainable  is  again  limited 
by  the  minimum  jet  impulse  available  from  the  thrusters. 
Figure  5-7  utilizes  the  third  switching  curve  and  can  be 
compared  to  Figure  5-6  where  only  two  switching  curves  were 
realized  to  demonstrate  the  improved  fuel  consumption. 


Figure  5-8  shows  the  time  responses  of  thruster  operation 
during  fine  tuning  in  both  cases.  Note  the  single  impulse 
firing  that  occurs  at  the  edge  of  the  limit  boundaries  when 
the  third  switching  curve  is  implemented.  The  appendix 
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provides  a  description  of  the  Tutsim  program  that  simulates 
the  independent  control  law  about  a  single  axis. 
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Figure  5-8  Fine  Control 

C.  ACQUISITION  PHASE  CONTROL 

The  acquisition  phase  is  the  most  critical  aspect  of 
satellite  operation  after  a  successful  launch.  During 
acquisition,  the  satellite  may  be  tumbling.  Tumbling 
translates  to  high  error  rates  which  will  cause  coupling  about 
at  least  two  axes.  In  this  phase,  a  minimum  time  solution  is 
desired  to  cancel  the  excess  rates  quickly.  [6]  The  Euler 
moment  equations  of  motion  are 
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(5.17) 


IXWX  =  (ly  “  Iz)  WyWx  +  Tx 

IyWy  =  (IZ  “  IX)WXWZ  +  Ty 

I*"z  =  (I*  ”  Iy)"x"y  +  T* 

where  Tx/  Ty  and  Tz  are  the  control  torques;  the  product  of 
thrust  and  moment  arm.  If  a  single  axis  of  symmetry  is 
assumed,  as  is  the  case  for  a  cylinder,  two  of  the  moments  of 
inertia  will  be  equal 

I*  -  I*  -  I 


The  moment  equations  become 

wx  =  a«yWz  +  Ux  (5.18) 

"y  “  Uy 

wz  —  -aUyWx  +  Uz 

where  a  =  (I-Iy)/I  and  u^Ti/Ii. 

Motion  about  the  axis  of  symmetry  is  independent  of  the 
other  two  axes  and  is  controlled  at  all  times  as  discussed  in 
the  previous  section.  Assuming  that  wy  is  constant  when  u,  is 
zero  during  rate  reduction,  the  state  equations  become 


'Xl(t)' 

'0  R 

Xl(t)‘ 

1  o' 

,X2(t)_ 

-R  0. 

,x2(t). 

0  1 

u(t) 


(5.19) 
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where  A =a«y/  xx=wx,  x2=wz,  Uj=ux  and  u2=uz.  The  minimum  time 
cost  function  is 

J  =  St  dt  (5.20) 

The  Hamiltonian  for  the  coupled  system  using  both  controllers, 
ux  and  u2,  is 

H(x,u,p,t)  =  1  +  pxx2  -  p2x:  +  p^  +  p2u2  (5.21) 

for  A=l.  The  costate  equations  are 

Pi(t)  =  -  3H/3x1  =  p2(t)  (5.22) 

p2(t)  =  -  8 H/ 3 x2  =  -Pi(t) 

which  have  solutions 

p:(t)  =  c3  cos(t+6)  (5.23) 

p2(t)  =  c4  sin (t+6 ) 

To  minimize  the  Hamiltonian  on  an  optimal  trajectory,  the  form 
of  optimal  control  is 

ux(t)  =  -sign(Pi)  (5.24) 

u2(t)  =  -sign  (p2) 
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The  control  torque  of  each  set  of  thrusters  will  reverse  every 
half  cycle  while  the  control  mix  will  change  every  quarter 
cycle.  The  minimum  time  switching  curves  can  be  derived  using 


dxx  _  x2  +  ux 
dx2  -Xi  +  u2 

/(-Xi+u^dXi  =  /(Xa+u^dXz 

2?L  “  u2xi  =  “2EL  ~  UjX2 
2  2 


(5.25) 


Adding  u22+Ui2  to  both  sides  of  (5.25)  and  multiplying  by  two, 
yields 


(Xj  -  u2)2  +  ( x2  +  uj2  =  ux2  +  u22  (5.26) 

For  normalized  control  torque,  u:=±l,  u2=±l,  (5.26)  becomes 

(x,  T  l)2  +  (x2  ±  l)2  =  2  (5.27) 

Equation  (5.26)  describes  four  circles  with  their  centers  at 
(±1,±1)  in  each  quadrant.  The  switching  curves  are  defined 
as  the  quarter  of  each  circle  that  intersects  the  abscissa  and 
ordinate  in  each  quadrant  as  shown  in  Figure  5-9.  Beyond 
these  scalloped  sections  the  switching  curves  are  mirror 
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images  of  the  half  sinusoidal  sections  along  each  axis,  but 
minimal  error  is  experienced  if  switching  is  implemented  about 
the  axes  themselves. 


l*x 

Figure  5-9  Double  Controller 
Minimum  Time  Switching  Curves 

A  phase  plane  representation  is  used  to  illustrate  rate 
reduction  where  wx  is  the  abscissa  and  wx  is  the  ordinate. 
Figure  5-10  depicts  the  implementation  of  the  four  quarter 
section  switching  curves  in  a  circular  region  with  a  radius 
of  two  for  initial  rate  conditions  of  (±1.5, ±1.5).  Figure  5- 
11  portrays  switching  about  the  coordinate  axes  in  the  region 
beyond  a  radius  two.  In  the  outer  region,  ux=-sign (wx)  and 
ux=-sign(wx)  .  Figure  5-12  shows  the  control  torque  time 
response  for  ux  and  ur. 
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Figure  5-12  Outer  Region  Control 


Since  the  system  is  coupled,  only  one  controller  is  needed 
to  reduce  the  rates  about  both  axes.  If  Uj=0,  equation  (5.25) 
becomes 

v  2  -v  2 

_i_  “  U2XX  =  _2  (5.28) 

2  2 


Adding  u22  to  both  sides  and  multiplying  by  two 

(*i  -  u2)2  +  x22  =  u22  (5.29) 

which  becomes 


( Xj  T  1 ) 2  +  X2Z  =  1 


(5.30) 
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for  u2=±l.  Conversely,  for  u2=0  equation  (5.30)  becomes 

Xj2  +  (X22  ±  1)  =  1  (5.31) 

for  ux=±l.  Equations  (5.30)  and  (5.31)  describe  two  symmetric 
circles  with  centers  at  x^fl  and  x2=±l,  respectively.  The 
switching  curves  for  the  minimum  time  solution  are  the  half 
circles  shown  in  Figure  5-13.  Figure  5-14  demonstrates  rate 
reduction  for  a  single  controller  system.  For  initial  rates 
in  the  outer  region,  switching  about  the  axes  was  implemented. 
Control  torque  values  of  u2=sign(ux)  for  (5.30)  and  ux=sign(w2) 
for  (5.31)  were  applied.  Figure  5.15  demonstrates  large  rate 
error  reduction  for  both  single  controller  cases. 


Curves 
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Figure  5-15  Outer 
Region  wT=wz=6 
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Figure  5-16  illustrates  the  control  torque  time  responses  m 
the  outer  region.  Comparison  of  Figures  5-12  with  5-16  shows 
that  the  double  controller  system  takes  half  the  time  using 
the  same  amount  of  control  effort  as  the  single  controller 
system. 


CONTROL  TIME  RESPONSE  Wx  =  W7.=6 
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Figure  5-16  Outer  Region  Control  Torque 


D.  LIMIT  CYCLE  CONTROL  FOR  ORION 

ORION  specifications  can  easily  be  implemented  in  the 
simulation  described  in  previous  sections.  A  three— axis 
stabilized  ORION  will  be  oriented  such  tnat  the  yaw  axis  is 
colinear  with  the  longitudinal  axis.  The  symmetry  of  the 
principal  moments  of  inertia  becomes  Ix=Iy  vice  Ix=Ir  as  in  the 
spin  stabilized  configuration.  ORION  principal  moments  of 


inertia  from  Table  4-1  are 


Ix  *  Iy  *  6.753  kg-m2 
Iz  =  2.776  kg-m2 

Attitude  control  thruster  data  is  obtained  from  Figure  4- 
14b.  The  thrust  available  from  the  attitude  control  reaction 
jets  is  between  .24  and  1.068  Newtons.  The  minimum  impuxse 
bit  is  .0223  N-sec  for  .025  seconds.  The  total  pulses 
available  as  a  function  of  the  catalyst  effective  lifetime  are 
161,000.  The  mass  flow  rate  of  the  thrusters  is  between  .104 
and  .453  grams/sec. 

Torque  is  the  product  of  thrust  and  moment  arm.  The 
length  and  radius  of  the  ORION  cylinder  are  .889  and  .2414 
meters,  respectively.  Thrusters  operate  in  pairs  to  prevent 
translational  movement  so  the  average  torque  produced  for 
ORION  is 


Tx  =  Ty  =  2(. 445)  (.889)  =  .7912  N-m 
Tz  =  2 ( .445) ( .2414)  =  .2148  N-m 

The  control  accelerations  are 

ux  =  My  =  .1172  s'2 
uz  =  .0774  s"2 
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Pointing  accuracies  of  ±1°  and  ±.05Vsec  were  used  in  the 
example  illustrated  to  define  the  limit  region.  The  minimum 
fuel-time  switching  curves  were  implemented  with  the  values 
of  A  previously  used,  A=2.25  and  250  for  the  outer  region 
deadzone  and  limit  region  deadzone  switching  curves, 
respectively.  Figures  5-17  and  5-18  illustrate  the  minimum 
fuel-time  control  law  for  yaw  and  roll -pitch  when  ORION  values 
for  torque  and  moments  oi  inertia  are  used. 


Figure  5-19  shows  the  yaw  thruster  operation  for  initial 
conditions  of  angle=.25  rad  and  zero  rate. 
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Figure  5-18  Roll-Pitch  Limit  Cycle  for 
Limits  of  ±1°,  ± . 05 “/sec 


It  is  difficult  to  predict  fuel  usage  up  to  the 
operational  phase  of  flight  for  three-axis  control.  The 
approximate  limit  cycle  lifetime  with  all  thrusters  operating 
is 


Total  Fuel  Available 

L  =  -  At  (5.32) 

Fuel  per  Control  Cycle 

where  At  =  one  limit  cycle.  [1]  In  the  limit  cycle,  with  no 
environmental  disturbances  added,  individual  pairs  of 
thrusters  pulse  on  four  times  per  cycle.  Using  a  mass  flow 
rate  of  .104  grams/s  for  .025  sec,  the  minimum  pulse  time, 
gives  a  total  fuel  per  control  cycle  of  between  .0624  and 
.0832  grams  if  all  six  thrusters  are  operating  to  maintain 
attitude  in  their  respective  limit  regions.  The  yaw  and  roll- 
pitch  control  cycle  duration  for  this  example  are  108  and  72 
seconds,  respectively.  Assuming  that  2/3  of  the  original 
32.35  kg  of  hydrazine  is  available  for  limit  cycle  attitude 
control,  total  limit  cycle  lifetime  would  be  about  320  days. 
Approximate  limit  cycle  lifetimes  for  various  combinations  of 
weighting  factors,  rate  limits  and  minimum  impulse  bits  based 
on  the  hydrazine  budget  are  listed  in  Table  5-1.  Increasing 
the  value  of  A  or  decreasing  the  rate  error  limit  will 
decrease  the  rate  error  and  thus  lengthen  the  thruster  off 
time,  but  the  maximum  effective  value  of  A  is  limited  by  the 
the  minimum  impulse  bit  on-period.  For  the  thruster  used 
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Table  5-1  LIMIT  CYCLE  LIFETIMES 


uz  Atz(sec)  rate(Vsec)  x  Atez(sec)  L(days) 


.0774 

.025 

BB 

48 

143 

.0774 

.01 

mm  m 

108 

323 

.0774 

.01 

HESS 

1000 

125 

427 

.0774 

.006 

.015 

1000 

185 

632 

.0387 

.025 

.075 

100 

86 

257 

.0387 

.025 

.075 

1000 

94 

281 

.0387 

.01 

.015 

1000 

218 

869 

u*y  Atxy(sec) 

rate ( '/sec) 

*  Atcxy  (sec) 

.1172 

.025 

.085 

30 

.1172 

.025 

.05 

72 

.1172 

.  007 

.025 

1000 

112 

.1172 

.004 

.015 

1000 

180 

.0586 

.025 

.075 

100 

58 

.0586 

.025 

.075 

1000 

64 

.0586 

.008 

.015 

1000 

195 

in  this  example,  A>100  has  negligible  effect  since  the  minimum 
AV  imparted  is 


AVZ  =  uz  At  =  .0774 (.025)  =  .001935  rad/sec 
AV^  =  Uxy  At  =  ,1172  (.025)  =  .00293  rad/sec 

which  will  kick  the  rate  error  from  the  negative  limit  of 
-8.73x10’*  rad/sec  to  a  positive  value  of  .0012  rad/ sec.  The 
example  illustrated  in  the  figures  uses  an  integration  step 
size  of  .01  seconds.  In  order  to  implement  the  fuel  saving 
choj  jes  listed  in  table  5-1  attitude  control  thrusters  with 
smaller  minimum  impulse  bits  and\or  on-times  must  be  utilized. 


Limit  cycle  lifetime  based  on  the  number  of  available 
pulses  ranges  from  30  to  90  days.  This  is  considerably  less 
than  limit  cycle  lifetimes  based  on  fuel  consumption.  In 
addition  to  smaller  impulse  bits,  the  number  of  available 
pulses  must  increase  to  lengthen  the  active  life  of  ORION. 
Another  option  would  be  to  operate  for  specified  periods 
during  each  orbit  vice  continuously.  Rather  than  replacing 
the  thrusters  with  more  complex  and  expensive  items,  this 
could  increase  the  active  lifetime  of  the  existing  design  to 
acceptable  levels. 

E.  ENVIRONMENTAL  DISTURBANCE  EFFECTS 

Since  the  majority  of  the  satellite  lifetime  will  be 
spent  in  the  operational  or  steady  state  mode,  this  is  the 
phase  when  the  effects  of  the  external  disturbances  of  gravity 
gradient  and  aerodynamic  drag  will  play  the  largest  part. 
Substituting  the  small  angle  approximation  for  the  angular 
velocity  vector  and  adding  a  disturbance  moment  to  the  Euler 
Moment  Equations  yields 


4>  =  a  w02  <f>  +  ux  +  Dx  (5.33) 

6  =  Uy  +  Dy 
=  uz  +  Dz 
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2  \  Iy"1* 

<t>  -  -2 u0  I  —  4> 

e  -  3w02  | 


Iz-I, 


(5.35) 


which  are  equations  of  motion  for  simple  harmonic  oscillators. 
The  moment  of  inertia  ratios  for  ORION  equal  +.5929  and 
-.5929,  respectively.  Figures  5-21  and  5-22  illustrate  the 
harmonic  motion  about  the  roll  and  pitch  axes.  Plots  in  5- 
21a  and  5-22a  are  phase  plane  representations  while  graphics 
in  5-21b  and  5-22b  show  the  time  response  for  angle  and  rate. 
The  motion  is  stable  and  periodic  with  a  period  that 
corresponds  to  the  satellite  orbit.  The  orbital  period  is 
defined  as  [4] 


P 


2  JT 


R3/2 


(5.36) 


Figure  5-23  furnishes  the  orbital  period  versus  altitude  for 
an  Earth  orbiting  satellite.  In  Figures  5-21  and  5-22  the 
maximum  rate  error  induced  never  exceeds  the  initial  value, 
but  the  angle  error  will  exceed  the  desired  pointing  accuracy 
over  the  course  of  an  orbit  if  a  rate  error  exists.  Gravity 
gradient  torque  has  a  time  constant  that  is  equivalent  to  an 
orbital  period.  The  time  constant  of  this  all-thruster 
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Figure  5-2lb  Gravitational  Roll  Effects 


Time  Response 
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Figure  5-22a  Gravitational  Pitch 
Effects  Phase  Plane 


Figure  5-22b  Gravitational  Pitch  Effects 
Time  Response 
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Figure  5-23  Orbital  Period  versus  Altitude 


attitude  control  system  is  on  the  order  of  100  seconds.  If 
the  control  system  is  operating  continuously  to  maintain  the 
desired  pointing  accuracy,  then  the  effects  of  gravity  are 
insignificant  compared  to  the  normal  cycling  that  occurs  in 
the  limit  region.  If  there  is  zero  rate  error  and  the  angle 
error  is  in  the  limit  region,  the  gravity  gradient  torque  will 
aid  in  maintaining  the  attitude  in  the  limit  region. 

The  disturbance  torques  contributed  by  aerodynamic  drag 

are 


L  Drag 


f  -Fd  L  sinflcos^  'j 
0 

Fd  L  COS0COS0 


(5.37) 
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where  FD  is  the  drag  force  as  a  function  of  atmospheric 
density.  ‘fable  5-2  lists  atmospheric  densities  at  some 
altitudes.  [8] 

Table  5-2  ATMOSPHERIC  DENSITIES 

Altitude (km)  Density (kg/m3) 


100 

6xl0'7 

200 

7xl0"7 

400 

lxlO'11 

1000 

lxlO"15 

At  h=400  km  the  drag  force  is  2.52xl0’<l  Newtons.  For  small 
angles  (5.37)  reduces  to 

tDx  =  ~Fd  L  8  =  -1 . 92x10  5  d  N-m 

TDz  =  Fd  L  =  1 . 92x10  5  N-m 

where  L=.076m  is  the  center  of  pressure  offset  from  the  center 
of  mass.  The  effects  of  aerodynamic  drag  alone  for  an  orbital 
period  are  illustrated  in  Figures  5-24  and  5-25  about  the  roll 
and  yaw  axes  for  initial  conditions  of  angle=.0l  rad  and 
rate=.0001  rad/sec.  The  pitch  angle  error  was  assumed 
constant  at  .01  rad.  The  satellite  response  to  drag  is 
unstable,  yet  like  gravity  gradient,  drag  causes  a  slow 
response  compared  to  the  continuously  operating  all-thruster 
control  system.  Drag  torques  that  would  cause  the  attitude 
of  the  satellite  to  be  unstable  over  the  course  of  an  orbital 
period  do  not  affect  motion  during  continuous  control. 
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Figure  5-25a  Drag  Yaw  Effects  Phase  Plane 


Finally  the  effects  of  gravity  gradient  and  aerodynamic 
drag  were  combined.  Substituting  the  disturbance  torques  into 
(5.33)  yields 


<t>  =  -a  u0z  <f>  +  \lx  -  2.84xl0'6  6  (5.38) 

which  is  the  equation  of  a  marginally  stable  oscillator.  The 
effect  that  gravity  and  drag  have  on  roll  is  illustrated  in 
Figures  5-26  and  5-27  for  initial  conditions  of  angle=.oi  rad, 
zero  rate  and  rate=.0001  rad/sec,  respectively.  The 
oscillation  does  not  increase  without  bound  and  actually  never 
repeats  the  same  cycle.  The  phase  plane  plots  have  a  run  time 
of  150  minutes  while  the  time  response  plots  have  a  run  time 
of  300  minutes.  Again,  the  time  constant  of  the  gravity-drag 
torque  is  orders  of  magnitude  bigger  than  that  of  the  all¬ 
thruster  control  system.  If  the  control  system  were  only 
operated  at  discrete  intervals  during  each  orbit  the  external 
disturbance  torques  could  together  run  the  rates  up  to 
unacceptable  levels. 
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Figure  5-26a  Gravity-Drag  Roll  Effects  with 
Zero  Rate  Error  Phase  Plane 


Figure  5-26b  Gravity-Drag  Roll  Effects  with 
Zero  Rate  Error  Time  Response 
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VI.  SUMMARY  AND  CONCLUSIONS 


The  purpose  of  this  thesis  was  to  design  an  all-thruster 
three-axis  attitude  control  system  for  ORION,  the  Naval 
Postgraduate  School  satellite  bus.  Pontryagin's  Minimization 
Principle  was  employed  in  the  derivation  of  fuel-time 
minimization  control  laws  for  the  acquisition  and  on-station 
phases  of  flight.  The  original  design  was  modified  slightly 
to  accommodate  a  three-axis  strategy. 

The  original  design  history  and  specifications  are 
discussed  in  Chapter  II.  The  total  mass  of  the  ORION 
satellite  is  250  lbm  including  71.5  lbm  of  the  monopropellant 
hydrazine  for  maneuvering.  The  attitude  control  system 
consists  primarily  of  one  5  lbf  thruster  for  orbital  transfer 
and  six  .1  lbf  thrusters  for  rotational  control.  The  original 
mass  placement  and  principal  moments  of  inertia  were  left 
unchanged  for  simulation  purposes.  The  four  stabilizing  booms 
are  not  utilized  in  this  analysis. 

Additional  concerns  that  affect  a  three-axis  stabilized 
satellite  compared  to  the  original  spin  stabilized  design  are 
addressed  in  Chapter  III.  To  meet  power  requirements,  solar 
panels  will  be  required  vice  just  body-mounted  solar  cells. 
A  three-axis  ORION  will  be  less  fuel  efficient  than  a  spinning 
version.  Various  schemes  using  different  combination*;  of 
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actuators  are  presented  for  future  study.  Alternate  booster 
selections  that  could  transfer  ORION  to  higher  orbits  than  the 
space  shuttle  and  minimize  the  use  of  the  40  lbf  thruster  are 
suggested . 

The  theory  behind  satellite  motion,  attitude 
determination  and  control  hardware  for  a  three-axis  stabilized 
spacecraft,  is  developed  in  Chapter  IV.  Rigid-body  spacecraft 
dynamics  and  the  effects  of  the  environmental  disturbances  of 
gravity  gradient  and  aerodynamic  drag  are  covered.  Attitude 
control  sensors  and  actuators  are  described  with  an  emphasis 
on  the  higher  degree  of  complexity  required  of  these 
instruments  vice  those  used  on  a  spinning  satellite.  Sensors 
include  horizon  sensors,  sun  sensors,  magnetometers,  and  rate 
gyros  while  the  actuators  are  strictly  thrusters. 

The  derivation  and  simulation  of  the  control  laws  that 
describe  the  acquisition  and  on-station  phases  of  flight  is 
found  in  Chapter  V.  it  was  determined  that  a  dual-axis 
minimum  time  control  law  was  optimal  for  reducing  the  rates 
as  quickly  as  possible  in  the  acquisition  phase.  Once  the 
rates  are  reduced  to  acceptable  levels,  a  fuel-time  control 
law  is  implemented  that  utilizes  a  weighting  factor  to 
minimize  fuel  usage. 

Two  regions  of  interest  developed  during  the  analysis  of 
the  fuel-time  control  law.  The  boundary  region  is  a 
rectangular  area  centered  about  the  desired  position.  The 
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deadzone  is  determined  by  the  nominal  pointing  accuracy  (±1° 
ORION  design  goal)  in  angle  and  the  minimum  thruster  impulse 
bit  for  rate.  The  boundary  region  is  distinguished  by  a  four 
pulse  limit  cycle  that  repeats  itself  every  50-200  seconds  for 
ORION.  Assuming  that  2/3  of  the  total  hydrazine  budget  is 
available  for  limit  cycle  maintenance,  lifetimes  of  4-28 
months  were  obtained  with  all  thrusters  operating.  The 
limiting  factor  for  the  selected  thrusters  was  the  number  of 
pulses  available  before  the  Shell  405  catalyst  lost 
effectiveness.  The  limit  cycle  lifetime  decreased  to  30-90 
days.  The  second  region  of  interest  was  outside  the  boundary 
region,  but  close  enough  for  the  satellite  to  still  experience 
independent  or  uncoupled  motion  about  one  or  more  axes.  A 
double  switching  curve  control  law  distinguished  by  "bang-off- 
bang"  control  was  realized  successfully. 

The  environmental  disturbances  of  gravity  gradient  and 
aerodynamic  drag  were  also  simulated  in  Chapter  V.  The 
gravity  gradient  disturbance  imitates  a  simple  harmonic 
oscillator  with  a  period  corresponding  to  an  orbital  period. 
A  typical  orbital  period  for  a  low  earth  orbit  is  on  the  order 
of  90  minutes.  Aerodynamic  drag  applies  an  unstabe  torque, 
but  when  coupled  with  gravity  causes  a  random  but  marginally 
stable  type  of  oscillation.  Both  disturbance  torques  are  so 
small  over  the  course  of  a  limit  control  cycle,  that  the 
effects  are  negligible. 
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The  three-axis  stabilized  design  simulated  in  this  thesis 
could  be  a  viable  attitude  control  option  for  the  general 
purpose  ORION  bus.  Upgraded  thrusters  must  be  selected  that 
have  smaller  minimum  impulse  bit  on-times  and  thrust  values 
in  addition  to  an  increased  catalyst  lifetime  if  ORION  is  to 
operate  in  the  constant  update  control  mode  for  a  feasible 
duration.  Another  possibility  that  is  only  mentioned  in  this 
work,  concerns  intermittant  operation.  ORION  could  function 
for  specified  periods  every  orbit  when  the  desired  target  is 
in  the  field  of  view.  This  method  would  insure  an  extended 
lifetime,  but  the  influence  of  external  torques  could  be  a 
limiting  factor. 
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APPENDIX 


ATTITUDE  CONTROL  AND  EXTERNAL  DISTURBANCE  PROGRAMS 


The  following  programs  simulate  a  three-axis  stabilized 
satellite  in  the  acquisition  and  on-station  phases.  The 
effects  of  gravity  gradient  and  aerodynamic  drag  are  modeled 
in  separate  routines. 

The  software  package  TUTSIM™  was  used  to  execute  the 
control  laws  derived  in  Chapter  V  [11].  TUTSIM™  models 
continuous  time  control  systems  in  a  signal  flow  format.  The 
following  commands  are  defined  for  the  potential  user: 


INT  -  integrator  block 
CON  -  constant  block 
GAI  -  gain  block 
MUL  -  multiplier  block 
DIV  -  divide  block 
ABS  -  absolute  value  block 
AND  -  logical  and  block 
OR  -  logical  or  block 
IFE  -  If-then  statement 
REL  -  relay  block 


The  numbers  preceding  commands  are  simply  block  identifiers. 
The  numbers  following  commands  are  block  inputs.  The  constant 
block  does  not  have  an  input.  The  integrator  and  gain  blocks 
require  at  least  one  input.  The  multiply,  divide  and  logical 
operator  blocks  need  at  least  two  inputs.  The  if-then 
statement  has  three  inputs.  The  first  input  is  the  argument. 
The  second  and  third  inputs  are  the  commands  that  support  true 
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or  false  arguments,  respectively.  The  relay  block  has  four 
inputs.  The  last  input  is  the  argument.  The  first  three 
inputs  are  commands  that  support  greater  than,  equal  to  or 
less  than  results,  respectively.  Further  information  on  the 
capabilities  of  TUTSIM™  can  be  found  in  the  operator  manual. 

The  first  program  simulates  on-station  control  about  one 
axis.  Input  parameters  include  initial  conditions,  principal 
moments  of  inertia,  control  torques,  limit  region  boundary 
values  and  weighting  factors  for  two  switching  curves.  The 
command  logic  applies  control  according  to  the  satellite 
attitude  in  the  phase  plane  in  relation  to  the  switching 
curves  and  boundary  limits.  The  control  law  for  this  program 
is  derived  in  the  first  section  of  Chapter  V. 

The  next  three  programs  simulate  the  acquisition  phase 
for  single  and  double  axis  control.  Input  parameters  include 
initial  conditions  for  rate  errors,  principal  moments  of 
inertia  and  torque  values.  These  programs  apply  control  as 
a  function  of  the  rates  for  coupled  equations  of  motion. 
Switching  occurs  about  the  ordinate  and  abscissa  in  a  defined 
outer  region.  Switching  occurs  on  minimum  time  curves  in  the 
inner  region.  The  control  law  for  the  acquisition  routines 
can  be  found  in  the  second  section  of  Chapter  V. 

The  last  two  programs  simulate  the  effects  of  gravity 
gradient  and  aerodynamic  drag  distubances,  respectively. 
Inputs  parmeters  are  provided  for  initial  conditions. 
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principal  moments  of  inertia  and  orbital  altitude.  The  last 
section  of  Chapter  V  illustrates  the  effects  of  external 
disturbances . 
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2  Acquistion  Control  for  a  Single 
Actuator  Set,  Uz 
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3  Acquisition  Control  for  a  Single 
Actuator  Set,  Ux 
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5  Gravity  Gradient  Disturbance  Effects 
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